Quantum Effects of Strain Influence on the Doping Energy in Semiconductors 
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Applying external strain is an efficient way to manipulate the site preference of dopants in semiconductors, 
however, the validity of the previous continuum elastic model for the strain influence on the doping forma- 
tion energy is still under debate. In this paper, by combining quantum mechanical theoretical analysis and 
first-principles calculations, we show that if the occupation change of different orbitals caused by the strain is 
negligible, the continuum elastic model is valid, otherwise it will fail. Our theory is confirmed by first-principles 
calculation of Mn-doped GaAs system. Moreover, we show that under compressive strain the hole density, thus 
the Curie temperature Tc can increase in Mn-doped spintronic materials. 

PACS numbers: 61.72.uj, 71.55.Eq, 71.15.Nc 
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Doping is an important technique to modify and tune mate- 
rial properties in semiconductors. (Ill Usually, efficient doping 
requires that the dopant stays at a specific site in the host crys- 
tal. For example, to improve the Curie temperature for the 
spintronics application of Mn-doped GaAs|2, 3], Mn atoms 
should substitute at the Ga site (Mnco) to create holes, and the 
existence of Mn at the interstitial site (Mn,) should be avoided 
because it acts as donors which will kill holes that mediate the 
ferromagnetism.lQ-Iil The site preference of a dopant is de- 
termined by its formation energies at different doping sites. 
Therefore, to increase Mnca and reduce Mn, concentration 
in GaAs, one should reduced the impurity formation energy 
of Mncfl and increase the impurity formation energy of Mn,. 
To achieve this, besides growing the sample at, e.g., Ga-poor 
condition, another approach could be applying strain (stress), 
because the change of dopant formation energy under strain is 
site dependent. [Wlw] 

Many studies and progresses have been made to understand 
how dopant formation energy at different sites is affected by 
strain. Zhang et al. have shown that an applied pressure can 
stabilize the substitutional position of transition metal impuri- 
ties in Si. 1 17] Theoretical calculations by Chen et al. indicate 
that the formation energy of self-interstitial atoms depends on 
strain. B 180 More recently, it has been shown by Zhu et al. us- 
ing first-principles calculation that, for an external hydrostatic 
or epitaxial strain, the change of the impurity formation en- 
ergy is monotonic, and it decreases if the external hydrostatic 
strain is applied in the same direction as the volume change 
caused by doping. il 9| In all these studies, classical continuum 
elastic theory iOl ll llll9l I20I1 is used to explain experimental 
and theoretical observations. In this theory, it assumes that the 
dopant induces a stress on the host crystal. When an external 
strain e is applied, the formation energy of the dopant, t^E, as 
a function of strain can be written as 



(1) 



where AEq is the dopant formation energy without strain, Vq 
is the equilibrium volume of the host. Pq is the pressure ex- 
erted on the host by the dopant. A similar relation between AE 
and external stress cr has also been proposed.! 9, 20] Eq. (1) 
clearly indicates that if the strain e has the opposite sign as 
Po, the dopant formation energy will decrease linearly as a 
function of e. The change will be large if V^Pd is large. 

This simple formula has been widely used in the past lfTll 
[T3UT5I1 . However, the origin of this formula is not well un- 
derstood, especially the validity of the linear relationship be- 
tween the strain and the formation energv. lfl2l [Tsl I2TI |22| Due 
to the importance of this formula for understanding the elastic 
theory, the site preference of dopants under strain, and the dif- 
fusion behavior between different sites, a detailed study and 
understanding are needed. 

In this paper, we find through quantum mechanical analy- 
sis and first-principles calculation that the precondition for the 
continuum elastic model and thus the linear dependence of the 
doping formation energy on strain is that there is no discon- 
tinuity in the occupation of the electron state when strain is 
applied. If the occupation of the energy level is switched after 
applying the strain due to band crossing in a quantum mechan- 
ical description, simple elastic model described by Eq. (1) will 
break down. This analysis is supported by our direct calcula- 
tion of the formation energy of Mn-doped GaAs at different 
doping sites. We find that the change of the formation energy 
of Mn, as a function of strain cannot fit to a linear equation 
because under compressive strain electrons from conduction 
band derived donor level transfer to the exposed unoccupied 
Mn d level. Moreover, we show that a compressive strain can 
enhance the concentration of free holes produced by M«g„ 
and suppress the formation of Mn,, thus increase the hole den- 
sity and consequently the Curie temperature Tc- 

We start our analysis by expanding the total energy of the 
host as a function of strain in a Taylor series. According to the 
quantum theory, for a solid under small strain (e), the energy 
up to the second order can be written as. 
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where Eq is the ground state energy of the host at e = 0, i/ and 
*P are the Hamihonian and wavefunction of the system, re- 
spectively, and the derivatives are computed at e = ^^y^ - 0, 
i.e., at y = Vq. Here, we assume volumetric strain. Similar 
expression for other type of strain such as epitaxial or uniaxial 
strain can also be derived in a similar way. The dopant forma- 
tion energy is the energy difference between the doped system 
and the undoped system. The total energy of the doped system 
can also be expanded into a formula like Eq. (1) except that Eq 
is replaced by Ed, the total energy of the doped system, and 
e' - ^y^" , the strain with respect to the equilibrium volume of 
the doped system. By expanding the series of the total energy 
of the doped system in terms of e instead of e' and subtracting 
that of the host total energy, we see that 



AE(e) = A^o + VoPd£ + ae 



(3) 



where, AEo - Eo - Eq is the dopant formation energy at equi- 
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v=Vd is the pressure exerted on the 



host by the dopant. The coefficient of the linear term VqPd 
is proportional to the volume difference 5Vd - Vd - Vq- It 
originates from the size difference between the dopant and the 
host elements. Because 5Vd is different for dopants at differ- 
ent sites, the slope of the change of AEi e) will be different for 
different dopants, thus this property can be used to tune site 
preference of a dopant. The coefficient a of the quadratic term 
depends on the elastic constant difference between the doped 
and undoped system, therefore, it is negligible if the doping 
concentration is small, i.e., Eq. (3) becomes Eq. (1) at the im- 
purity doping limit and if the applied strain is small. This 
explains the origin of the continuum elastic model of Eq. (1). 

However, in our derivation above, we assume that the 
change of the Hamihonian and wavefunction is a continuous 
function of e. This assumption is usually valid for semicon- 
ductors because the valence band and conduction band will 
not cross each other under strain before the metal-insulator 
transition. This situation also holds in most conventional 
doped case when the defect levels (e.g., Sica donor level in 
GaAs) is derived from the band edge state, so no band cross- 
ing is expected under strain. However, this assumption will 
not hold in a quantum systems where electrons can jump from 
one level to another level with very different atomic wavefunc- 
tion characters when strain is applied, thus cause a discontin- 
uous change in the Hamihonian, i.e., the second term of Eq. 
2. This usually happens when the defect level has different 
origin from the host energy level, e.g., the d orbital levels in 
a transition metal doped conventional semiconductors. From 
the analysis above, we can expect that after the charge trans- 
fer, the dopant changed from one state to another state with 
different size, therefore, 5Vd will change to 6Vd*- In this case, 
VqPd itself is also dependent on e, thus the linear dependence 
will break down. If the charge transfer is gradual, then there 
will be a gradual transition region where the slope changes. 

To demonstrate the above theory, we have calculated the 
formation energy of Mn, doping in GaAs at two interstitial 
sites and MnQa as a function of volumetric strain. Mn 3d band 
is located near the band edge of GaAs, which provides a pro- 
totype example to verify our theory above. The Mn-doped 
GaAs system is also a typical diluted magnetic semiconductor 



for the development of the spintronic devices, thus, the calcu- 
lated results also has significant technical importance^Iltl- 

The first-principles calculations are carried out using Vi- 
enna ab initio simulation program (VASP),|23| based on den- 
sity functional theory with the projector augmented wave 



(PAW) pseudopotential |24 
approximation (PBE-GGA) 



and the generalized gradient 
1 2511 to the exchange-correlation 
functional. The energy cut-off is set as 400 eV. The total en- 
ergy and force on the atoms are converged within 10 eV 
and 10"^ eV/A, respectively. A reciprocal space k-point mesh 
of 4x4x4 for all the supercells is employed. A (2x2x2) cu- 
bic supercell containing 64 host atoms is adopted for all the 
defect calculations. The optimized lattice constant is 5.75 A 
for the zinc -blende structure of pure GaAs, which is in good 
agreement with previous studies. ll26ll27ll 

In Fig. 1 we plot the strain influence on the calculated for- 
mation energy of Mn doping in GaAs at four different sites: 
Mni_As interstitial surrounded by four As atoms, Mriijja inter- 
stitial surrounded by four Ga atoms, Mn substitution at the Ga 
site, Mrica, and a defect complex 2Mnca + Mni r;/,.ll6l l28ll29ll 
The formation energy for a dopant at site a with the charge 
state q is defined as 



AH(a, q) = E{a, q) - E{GaAs) + nca[E(Ga) + ;u(Ga)] 
-nM« [^(Mn) + iJ(Mn)] + q[evBM + Ef], 



(4) 



where E(a, q) is the calculated total energy of the defect » in a 
supercell in the charge state q, EiGaAs) is the energy of pure 
GaAs in the same supercell. n, is the number of Ga and Mn 
atoms that are removed from the system during doping {nca - 
and UMn - -1 for the interstitial doping and nca - 1 and 
"Mn = -1 for the substitutional defect). ^(Ga) and ^(Mn) are 
the energy per atom of bulk Ga or Mn respectively, and ju(Ga) 
and yu(Mn) are chemical potential of Ga and Mn, referenced 
to the energy of bulk Ga and Mn, respectively. To improve p- 
type doping, Ga poor growth condition with ju(Ga)=-0.74 eV 
and /i(Mn)=-0.61 eV are used in our calculation. [28| evBM is 
the eigenenergy of the valence band maximum (VBM) state 
of pure GaAs, and Ef is the Fermi energy level referenced to 
the VBM energy. 

From Fig. 1, we can see two obvious differences in the 
formation enthalpy change of substitution and interstitial de- 
fects, (i) the formation energy increases with strain for Mnca 
, whereas decreases for Mtiij^s, Mni_Ga and the defect com- 
plex. This can be explained according to the volume change 
caused by the impurity formation. We find from our calcula- 
tions that the volume is decreased when a Mnca substitution 
is formed, and increased when the interstitials are formed, 
so increasing e will decrease the formation energy of inter- 
stitials and increase that of substitutional dopant. |19| (ii) the 
change is almost linear for Mtica and the defect complex from 
-0.1 15 < e < 0.125, in good agreement with that predicated 
by Eq. (1), whereas the change is not linear for the two in- 
terstitial cases, Mrtijis and Mni_ca, and is not agreement with 
what predicated by Eq. (1). In the following, we will explain 
why the nonlinearity exist for the two donor defects Mn,_A.s 
and Mtiijrja from the band structure of Mn-doped GaAs. 

As a typical semiconductor, it is known that the top of va- 
lence band of GaAs has mainly the bonding component of 



3 



5.0 B----I;---.. 

4.5 - 



S __, I I I I ■ 

111 ' ' ' ' ' 

£ i °r::: = * \ 

0.5 - 

-0.10 -0.05 0.00 0.05 0.10 

Volumetric Strain (e) 

FIG. 1. (Color online) The calculated formation energies as functions 
of the strain 6 for Mn doping in GaAs at different sites: Mn,_^s and 
Mrii Oa interstitials, Mnca substitution, and IMnga + Mtii Ca complex 
per Mn dopant. The predictions calculated from Eq. (1) (straight 
lines) are plotted for comparison. 




FIG. 2. (Color online) Schematic plot of the band components of 
the Mn-doped GaAs and the position of the Fermi energies of MriQa 
(left) and Mf!,_^j (right). 

the p-p hybridization between As and Ga, while the bottom 
of the conduction band mainly has the antibonding compo- 
nent of the s-s hybridization between As and Ga. When Mn, 
which has d^s^ atomic configuration, is introduced into GaAs, 
the spin-up Mn 3d states are fully occupied whereas the spin- 
down Mn 3d states are empty, producing impurity levels near 
the conduction band edges. 

For MnQa acceptor, it creates a hole at the top valence band, 
so without strain Mnca has cf^ + h configuration (Fig.2(left)). 
Because both VBM and the localized Mn d orbital has low 
absolute deformation potential, 130] when strain is applied, the 
order of the occupied energy level does not change. This ex- 
plains why the linear relationship between the formation en- 
ergy and strain is kept for Mnca, and can be well predicated 
by Eq. 1. 

In contrast to Mnca, Mrii acts as a donor, generating two 
electrons at the donor level near the conduction band max- 
imum (CRM). The spin-down Mn-3d band is only slightly 
higher than the donor level for Mn,-, so it is not occupied when 
there is no strain. Therefore, at e = 0, M«, has d^^ +2e config- 
uration (Fig.2(right)). When a compressive volumetric strain 



is applied to the system, the antibonding conduction band as 
well as the derived donor level with As 4s and Ga 4s char- 
acter shifts upward in energy. On the contrary, Mn 3d states 
are quite localized, so its energy level is only weakly influ- 
enced by the strain. Due to the upward shift of the donor level 
relative to the Mn-3d spin-down band, when e < -0.058 elec- 
trons start to transfer from the CBM derived donor level to the 
Mn-3d spin-down band, so the Mn,- configuration changed to 
d^^d^^ as the compressive strain is applied to the M«, doped 
system. The change of the atomic wavefunction character 
from s, p to d make the simple linear relation of Eq. (1) fail. 
Fig. 1 shows that under compressive strain, a linear line for 
e > is replaced by another linear line for e < -0.058 with 
a transition region between -0.058 < e < 0.0. The smaller 
slope associated with the high compressive strain is because 
Mn in d^^d^^ configuration at high compressive strain is more 
localized, i.e., has a smaller size, than in the d^^ + 2e configu- 
ration at zero or expansive strain. Similarly, the smaller slope 
for Mnijja than for Mnij^s is because there are more electrons 
around the anion atom As, so the pressure exerted by Mn at 
an interstitial site next to As is larger than that at an intersti- 
tial site next to Ga. All these are consistent with our theory 
discussed above. 

Our discussion above show that the reason that the for- 
mation energy of Mn, does not follow the simple linear re- 
lation of Eq. (1) is because Mn, has two electrons near the 
conduction band edge and the strain changes the electronic 
occupation of different bands. Therefore, if we form the 
Mriijja + 2Mnca complex, so the two electrons from Miiijja 
passivate the holes of Mrtca, then the strain in the range - 
0.115 to 0.125 will not change the electronic occupation of 
the bands, thus the dependence of its formation energy on the 
strain should be more linear Indeed our calculations con- 
firmed this expectation, as shown in Fig. 1. Therefore, our 
calculations demonstrated that in Mn-doped GaAs system, 
whether the doping formation energy changes linearly with 
the strain, i.e., the validity of the continuum elastic model, 
depends on whether the strain changes the electronic occupa- 
tion or not. It should be mentioned that, although the current 
analysis are demonstrated for the formation energy, the anal- 
ysis also apply to the energy differences and diffusion bamers 
between different doping sites. 

The above discussed dependence of formation energy on 
strain can be used to tune doping properties for a specific ap- 
plications. For example, the experimentally observed feiTO- 
magnetism of Mn-doped GaAs is mediated by the hole pro- 
duced by the Mtica acceptors. |6- 8] However, as the hole den- 
sity increases and the Fermi energy shifts towards VBM, com- 
pensating Mn, starts to form, which will lower hole density 
and thus the Curie temperature Tc- Our calculated results in 
Fig. 1 suggest that the relative ratio between Mnca and Mn,- 
will increase significantly when compressive strain is applied, 
therefore, the Tc should increase if compressive strain can be 
applied during the doping process. 

In conclusion, the strain influence on the doping forma- 
tion energy at different sites in semiconductors is analyzed 
in terms of the quantum mechanics theory and demonstrated 
by performing first-principles calculation of Mn-doped GaAs 
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system. We show that, the validity of the continuum elastic 
model, i.e., whether the formation energy is a linear function 
of the applied strain, depends on whether the occupation of 
the difference electronic bands is changed or not. If the oc- 
cupation change caused by the strain is negligible, the forma- 
tion energy is linearly dependent on the strain, while if the 
occupation change is significant, the linear relation, thus the 
continuum elastic model will fail. Our study clarifies the pre- 
vious confusion about this linear relationship, and provides an 
easy way for predicting the strain influence on the doping site 



in semiconductors. The calculation also shows that the Fermi 
energy pinning level of Mn-doped GaAs can be shift to close 
to VBM through applying compressive strain to the system, 
which can increases the ferromagnetic Curie temperature. 
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